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Abstract 

The Degasperis-Procesi equation can be derived as a member of a one- 
parameter family of asymptotic shallow water approximations to the Eu- 
ler equations with the same asymptotic accuracy as that of the Camassa- 
Holm equation. In this paper, we study the orbital stability problem of the 
peaked solitons to the Degasperis-Procesi equation on the line. By con- 
structing a Liapunov function, we prove that the shapes of these peakon 
solitons are stable under small perturbations. 
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1 Introduction 

The Degasperis-Procesi (DP) equation 

y t + y x u + 3yu x = 0, x e R, t > 0, (1.1) 

with y = u — u xx , was originally derived by Degasperis-Procesi [H] using the 
method of asymptotic integrability up to third order as one of three equations 
in the family of third order dispersive PDE conservation laws of the form 

u t - a 2 u xxt + ju xxx + c u x = (ciu 2 + c 2 u x + c 3 uu xx ) x . (1.2) 

The other two integrable equations in the family, after rescaling and applying a 
Galilean transformation, are the Korteweg-de Vries (KdV) equation 

u t + u xxx + uu x = 

and the Camassa-Holm (CH) shallow water equation [2l [Pol [TBI [22] . 

Vt + y x u + 2yu x = 0, y = u-u xx . (1.3) 
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These three cases exhaust in the completely integrable candidates for (|1.2[) by 
Painleve analysis. Degasperis, Holm and Hone [13] showed the formal integra- 
bility of the DP equation as Hamiltonian systems by constructing a Lax pair 
and a bi-Hamiltonian structure. 

The Camassa-Holm equation was first derived by Fokas and Fuchassteiner 
|18j as a bi-Hamiltonian system, and then as a model for shallow water waves 
by Camassa and Holm [2]. The DP equation is also in dimensionless space- 
time variables (x, t) an approximation to the incompressible Euler equations 
for shallow water under the Kodama transformation [H] and its asymptotic 
accuracy is the same as that of the Camassa-Holm (CH) shallow water equation, 
where u(t, x) is considered as the fluid velocity at time t in the spatial x-direction 
with momentum density y. 

Recently, Liu and Yin [24] proved that the first blow-up in finite time to 
equation (II. ip must occur as wave breaking and shock waves possibly appear 
afterwards. It is shown in [24] that the lifespan of solutions of the DP equation 
(jl.ip is not affected by the smoothness and size of the initial profiles, but affected 
by the shape of the initial profiles (for the CH equation, see [HIE])- This can be 
viewed as a significant difference between the DP equation (or the CH equation 
) and the KdV. It is also noted that the KdV equation, unlike the CH equation 
or DP equation, does not have wave breaking phenomena [30] . Under wave 
breaking we understand that development of singularities in finite time by which 
the wave remains bounded but its slope becomes unbounded |31j . 

It is well known that the KdV equation is an integrable Hamiltonian equation 
that possesses smooth solitons as traveling waves. In the KdV equation, the 
leading order asymptotic balance that confines the traveling wave solitons occurs 
between nonlinear steepening and linear dispersion. However, the nonlinear 
dispersion and nonlocal balance in the CH equation and the DP equation, even 
in the absence of linear dispersion, can still produce a confined solitary traveling 
waves 

u(t, x) = ap(x — ct) (1-4) 

traveling at constant speed c > 0, where <p(x) = e - '^'. Because of their shape 
(they are smooth except for a peak at their crest) , these solutions are called the 
peakons [2j [13] . Peakons of both equations are true solitons that interact via 
elastic collisions under the CH dynamics, or the DP dynamics, respectively. The 
peakons of the CH equation are orbitally stable [12] . For waves that approximate 
the peakons in a special way, a stability result was proved by a variation method 

Note that we can rewrite the DP equation as 

u t - u txx + 4:uu x = 3u x u xx + uu xxx , t > 0, x e M. (1.5) 

The peaked solitons are not classical solutions of (|1.5p . They satisfy the 
Degasperis-Procesi equation in the conservation law form 

u t + d x ():U 2 + )-ip* ( lu 2 ) J =0, t>0, (1.6) 
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where * stands for convolution with respect to the spatial variable This 
is the exact meaning in which the pcakons are solutions. 

Recently, Lundmark and Szmigielski [26) presented an inverse scattering ap- 
proach for computing n-peakon solutions to equation (|1.5|) . Holm and Staley [20) 
studied stability of solitons and peakons numerically to equation (|1.5p . Anal- 
ogous to the case of Camassa-Holm equation jBJ, Henry [115] and Mustafa [2T)| 
showed that smooth solutions to equation (|1.5p have infinite speed of propaga- 
tion. 

The following are three useful conservation laws of the Degasperis-Procesi 
equation. 

E\(u) = I ydx, Ei(u) = I yvdx, E%{u) — I u dx, 
Jr Jr Jr 

where y = (1 — d 2 )u and v — (4 — & x v )~ 1 u, while the corresponding three useful 
conservation laws of the Camassa-Holm equation are the following: 

Fi(u)= / ydx, F 2 (u)= / {u 2 +u 2 x )dx, F 3 (u) = / {u 3 +uu 2 x )dx. (1.7) 
Jr Jr Jr 

The stability of solitary waves is one of the fundamental qualitative proper- 
tics of the solutions of nonlinear wave equations. Numerical simulations [T31I25] 
suggest that the sizes and velocities of the peakons do not change as a result of 
collision so these patterns are expected to be stable. Furthermore, it is observed 
that the shape of the peakons remains approximately the same as time evolves. 
As far as we know, the case of stability of the peakons for the Camassa-Holm 
equation is well understood by now [TH [T2] , while the Degasperis-Procesi equa- 
tion case is the subject of this paper. The goal of this paper is to establish a 
stability result of peaked solitons for equation (|1.5[) . 

It is found that the corresponding conservation laws of the Degasperis- 
Procesi equation are much weaker than those of the Camassa-Holm equation. 
In particular, one can see that the conservation law ^(u) for the DP equation 
is equivalent to ||u||^2- In fact, by the Fourier transform, we have 

E 2 (u) = [ yvdx = / i±|J|^)| 2 ^ ~ INII* = (1-8) 

Therefore, the stability issue of the peaked solitons of the DP equation is more 
subtle . 

For the DP equation, we can only expect to obtain the orbital stability of 
peakons in the sense of L 2 — norm due to a weaker conservation law E^- The 
solutions of the DP equation usually tend to be oscillations which spread out 
spatially in a quite complicated way. In general, a small perturbation of a 
solitary wave can yield another one with a different speed and phase shift. We 
define the orbit of traveling-wave solutions cip to be the set U((f) = {cip(- + 
xo), xq GM.}, and a peaked soliton of the DP equation is called orbitally stable 
if a wave starting close to the peakon remains close to some translate of it at 
all later times. 
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Let us denote 

= h\\x- 

The following stability theorem is the principal result of the present paper. 

Theorem 1 (Stability) Let cip be the peaked soliton defined in \1.J$ traveling 
with speed c > 0. Then cip is orbitally stable in the following sense. If uq G H s 
for some s > 3/2, yo = uq — d^iio is a nonnegative Radon measure of finite total 
mass, and 

„ 1 

||"o - cip\\x < ce, \E 3 (u ) - E s (ap)\ < c } e, < e < -, 

then the corresponding solution u(t) of equation \1.5}) with initial value u(0) = uq 
satisfies 

sup\\u(t,')-ap(--Z 1 (t))\\ x <3c e 1/4 , 
t>o 

where £,i(t) G K is the maximum point of the function v(t, •) = (4 — u{t, •). 
Moreover, let 

Mi (t) = v(t, ^i(t)) > M 2 (*)•■■> M n (t) >Qandm 1 (t)>---> m„_i (*) > 

be all local maxima and minima of the nonnegative function v(t, •), respectively. 
Then 

Mi (t) - - < cV2e (1.9) 
6 

and 

n 

Y,(MHt)-rnl 1 (t))<2c^e. (1.10) 

Remark 1 The state of affairs about these maxima/minima implied by the pre- 
vious theorem is a consequence of the assumption on yo, as shown in Lemma 
\3.1[ For an initial profile uq £ H s , s > 3/2, there exists a local solution 
u G C([0, T),H S ) of \1.5}) with initial data u(0) = uo Under the assump- 

tion yo = uq — O^uq > in Theorem[l\ the existence is global in time Y<H\ , that 
is T = +oo. For peakons cip with c > 0, we have (l — <9^) (c(p) = 2c8 (here S 
is the Dirac distribution). Hence the assumption on yo that it is a nonnegative 
measure is quite natural for a small perturbation of the peakons. Existence of 
global weak solution in H 1 of the DP equation is also proved in [16]. Note that 
peakons cip are not strong solutions, since ip G H s , only for s < 3/2. 

The above theorem of orbital stability states that any solution starting close 
to peakons cip remains close to some translate of cip in the norm \\ \\x, oi 
any later time. More information about this stability is contained in il.9\) and 
il.lO\) . Notice that for peakons cip, the function v ctp is single-humped with the 
height ^c. So \1.9\) and hl.lCtjl imply that the graph ofv(t,-) is close to that of 
the peakon cip with a fixed c > for all times. 
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There are two standard methods to study stability issues of dispersive wave 
equations. One is the variational approach which constructs the solitary waves 
as energy minimizers under appropriate constraints, and the stability automat- 
ically follows. However, without uniqueness of the minimizer, one can only 
obtain the stability of the set of minima. The variational approach is used in 
PTj for the CH equation. It is shown in [TT] that the each peakon ctp is the 
unique minimum (ground state) of constrained energy, from which its orbital 
stability is proved for initial data uq £ H 3 with y = (1 — d 2 )u > 0. Their 
proof strongly relics on the fact that the conserved energy F2 in (jTTTJ) of the CH 
equation is the H 1 — norm of the solution. However, for the DP equation the 
energy E2 in (|1.8[) is only the L 2 norm of the solution. Consequently, it is more 
difficult to use such a variational approach for the DP equation. 

Another approach to study stability is to linearize the equation around the 
solitary waves, and it is commonly believed that nonlinear stability is governed 
by the linearized equation. However, for the CH and DP equations, the non- 
linearity plays the dominant role rather than being a higher-order correction to 
linear terms. Thus it is unclear how one can get nonlinear stability of peakons 
by studying the linearized problem. Morover, the peaked solitons ctp are not 
diffcrentiable, which makes it difficult to analyze the spectrum of the linearized 
operator around ctp. 

To establish the stability result for the DP equation, we extend the approach 
in [T^] for the CH equation. The idea in [T2] is to directly use the energy F2 as 
the Liapunov functional. By expanding F2 in (|1.7p around the peakon ctp, the 
error term is in the form of the difference of the maxima of ctp and the perturbed 
solution u. To estimate this difference, they establish two integral relations 



and the error estimate \M — max tp\ then follows from the structure of the above 
polynomial inequality. 

To extend the above approach to nonlinear stability of the DP peakons, we 
have to overcome several difficulties. By expanding the energy E2 (u) around 
the peakon ctp, the error term turns out to be maxv cip — maxu u , with v u = 
(4 — d 2 )~ 1 u. We can derive the following two integral relations for M\ = max?),,, 
E 2 (u) and E 3 (u) by 



with some functions g and h related to v u . To get the required polynomial 
inequality from the above two identities, we need to show h < 18maxw u . How- 
ever, since h is of the form —d 2 v u ± Qd x v u + 16w u , generally it can not be 




and 




with a function g. Relating these two integrals, one can get 



F 3 (u) < MF 2 (u) - -M 3 , M = maxu(x) 
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bounded by v u . This new difficulty is due to the more complicated nonlinear 
structure and weaker conservation laws of the DP equation. To overcome it, we 
introduce a new idea. By constructing g and h piecewise according to mono- 
tonicity of the function v u , we then establish two new integral identities (|3.7p 
and (|3.9p for £"2, £3 and all local maxima and minima of v u . The crucial es- 
timate h < 18maxw„ can now be shown by using this monotonicity structure 
and properties of the DP solutions. This results in inequality (|3.13[) related to 
£2, £3 and all local maxima and minima of v u . By analyzing the structure of 
equality (|3.13[) . we can obtain not only the error estimate \M\ — maxw cv | but 
more precise stability information from (|1.10[) . We note that the same approach 
can also be used for the CH equation to gain more stability information (see 
Remark [2) . 

Although the DP equation is similar to the CH equation in several aspects, 
we would like to point out that these two equations are truly different. One of 
the novel features of the DP equation is it has not only peaked solitons [13] . 
u(t,x) = ce~\ x ~ ct \ c > but also shock peakons [HHS] of the form 

u(t,x) = — sgn(x)e~\ x \, k > 0. (1-H) 

t + k 

It is noted that the above shock-peakon solutions [25] can be observed by sub- 
stituting [x, t) 1 — > (ex, et) to equation l|1.5[) and letting e — > so that it yields 
the "derivative Burgers equation" (u t + uu x ) xx — 0, from which shock waves 
form. The periodic shock waves were established by Escher, Liu and Yin [17]. 

The shock peakons can be also observed from the collision of the peakons 
(moving to the right) and antipeakons (moving to left) [2"5] . 

For example, if we choose the initial data 

uo(a:) = Cl e^ x - x ^ -cie-l"-**', 

with ci > 0, and x± + x-i — 0, X2 > 0, then the collision occurs at x — and the 
solution 

u(x,t) = Pl (t)e-^- gi ^ +p 2 (t)e^ x - q2< - t)l , 

(x,t) E R + xR, only satisfies the DP equation for t < T. The unique continuation 
of u(x, t) into an entropy weak solution is then given by the stationary decaying 
shock peakon 

On the other hand, the isospectral problem in the Lax pair for equation ()1.5|) 
is the third-order equation 

i> x - ^xxx - M)^ = 

cf. [T3j, while the isospectral problem for the Camassa-Holm equation is the 
second order equation 
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(in both cases y = u — u xx ) cf. [2]. Another indication of the fact that there is 
no simple transformation of equation (|1.5p into the Camassa-Holm equation is 
the entirely different form of conservation laws for these two equations [5J |T3] . 
Furthermore, the Camassa-Holm equation is a re-expression of geodesic flow on 
the diffeomorphism group [5J [TU] and on the Bott-Virasoro group [HI [M| , while 
no such geometric derivation of the Degasperis-Procesi equation is available. 

The remainder of the paper is organized as follows. In Section 2, we recall 
the local well-posedness of the Cauchy problem of equation (|1.5[) . the precise 
blow-up scenario of strong solutions, and several useful results which are crucial 
in the proof of stability theorem for equation (|1.5[) from [33J [33J . Section 3 is 
devoted to the proof of the stability result (Theorem [T]) . 

Notation. As above and henceforth, we denote by * convolution with respect to 
the spatial variable x G R. We use || • \\lp to denote the norm in the Lebesgue 
space L P (R) (1 < p < oo), and || • \\h=, s > for the norm in the Sobolev spaces 
H S (R). 



2 Preliminaries 

In the present section, we discuss the issue of well-posedness. The local exis- 
tence theory of the initial-value problem is necessary for our study of nonlinear 
stability. We briefly collect the needed results from [MJ [35J [33J . 

Denote p(x) := \e~\*\, x G R, then (1 - c^)" 1 / = p * f for all / G L 2 (R) 
and p * (u — u xx ) = u. Using this identity, we can rewrite the DP equation (| 1 . 5|) 
as follows: 

u t + uu x + d x p* (^u 2 ^ = 0, t > 0, x e R. (2.1) 

The local well-posedness of the Cauchy problem of equation (|1.5[) with initial 
data Mo G i/ 5 (R), s > | can be obtained by applying Kato's theorem [231 [32] . 
As a result, we have the following well-posedness result. 

Lemma 2.1 JffffJ / Given u G H S (M), s > |, t/iere exist a maximal T — 
T(uo) > and a unique solution u to equation hi. 5}) (or equation \2.1)) ). such 
that 

u = u(., uq) G C([0, T); iJ s (R)) n ^([O, T); H S ~ 1 (M.)). 

Moreover, the solution depends continuously on the initial data, i. e. the mapping 
u i ► u(-,uo) : -ff s (R) -> C*([0,r) ;j ff s (R)) n C 1 ([0,T);i? s - 1 (R)) is continuous 
and the maximal time of existence T > can 6e chosen to be independent of s. 

The following two lemmas show that the only way that a classical solution 
to (|1.5p may fail to exist for all time is that the wave may break. 

Lemma 2.2 f32f Given uo G H S (R), s > |, Wow up o/ t/ie solution u = u(-,Uo) 
in finite time T < +oo occurs if and only if 

liminfj inf \u x (t, x)]\ — — oo. 
tfT xeffi 
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Lemma 2.3 \2$ Assume uq £ H s (M.), s > |. Let T be the maximal existence 
time of the solution u to equation il.5\) . Then we have 



\\u{t,x)\\ L c- <3||uo(x)|||at+||tio(a!)|U», Vie [0,7]. 

Now consider the following differential equation 

q t = u(t,q), t £ [0,T), , s 

q(0,x) = x, x £ R. 1 ' > 

Applying classical results in the theory of ordinary differential equations, 
one can obtain the following two results on q which are crucial in the proof of 
global existence and blow-up solutions. 

Lemma 2.4 f33f Let uq £ H S (R), s > 3, and let T > be the maximal exis- 
tence time of the corresponding solution u to equation i 1.5]) . Then the equation 
12. lp has a unique solution q £ C ([0, 7) x R, R). Moreover, the map q(t, •) is 
an increasing diffeomorphism of R with 

q x (t, x) = exp (J u x (s, q{s, x))ds^j > 0, V(t, x) £ [0, T) X R. 

Lemma 2.5 \33f Let uo £ H S (M.), s > 3, and let T > be the maximal 
existence time of the corresponding solution u to equation H2.2\) . Setting y :— 
u — u xx , we have 

y{t,q(t,x))q 3 x (t,x)=y (x), V{t,x) £ [0,7) x R. 

The next two lemmas clearly show that the solution of equation (|1.5j) is 
affected by the shape of the initial profiles, not the smoothness and size of the 
initial profiles. 

Lemma 2.6 \21$ Let uq £ H S (M.), s > I. Assume there exists xq £ R such that 

y (x) = u Q (x) - u . xx (x) > if x<x , 

y (x) =u (x) -u , xx (x) < if x>x , 

and j/o changes sign. Then, the corresponding solution to equation \1.5\) blows 
up in a finite time. 

Lemma 2.7 \2$ Assume uq £ H S (M.), s > | and there exists xo £ R sucft i/iai 

y (x) < if x <x , 

y (x) > if x >x . 

Then equation il.5}) has a unique global strong solution 

u = u(.,«o) € C([0, oo); # S (R)) n C^QO, oo); H'^QS)). 

Moreover, E2(u) = J m yvdx is a conservation law, where y = (1 — <9^)u and 
v = (4 — and for all t £ R + we /iaue 

(%) ti*(V) > -KV)I onM > 

(ii) \\u\\i <6||uo||| 2 i 2 +4|| U n|li2||uolU~i+||uo||?- 
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The following lemma is a special case of Lemma 12.71 

Lemma 2.8 J£$ Assume u £ H S (R), s > §. If y = u - u , xx > (< 0) on 
R, then equation \1.5}) has a unique global strong solution u such that u(t,x) > 
(< 0) and y(t, x) = u - d 2 x u > (< 0) for all (t, x) £R + x E. 

Lemma 2.9 Assume uq £ H s , s > 3/2 and yo > 0. If k\ > 1, then the 
corresponding solution u of il.5\) with the initial data u(0) — uq satisfies 

(k 1 ±d x )u(t,x)>0, V(M)£l+xI. 

Proof. By Lemma [2.41 and a simple density argument, it suffices to show the 
lemma for s = 3. In view of Lemma [2.51 the potential y(t,x) = (1 — d 2 )u > 
0, \f(t,x) £ R + x R. Note u = (1 - d 2 )~V Then we have 

—x rx x poo 

u(t,x) = —J^ e r 'y(t,v)dv + Y J e~"«(*> »?)*7 (2.3) 

and 

e -x px e x poo 

u x (t,x) = —J e r 'y{t 1 i 1 )dr ! + — J e-*>y(t,r))dr). (2.4) 

It then follows from the above two relations (|2.3[) and (|2.4p that 



(fci ± d x )u = T l)e~ a / e"yd?7 + hk t ± l)e" / 



e-iydr) > 0. (2.5) 



□ 



Lemma 2.10 Let w(t,x) — (k\ ± d x )u(t,x). Assume uo <E -ff s , s > 3/2 and 
yo > 0. J/fci > 1 and fc 2 > 2, f/ien we /lave 

(fc 2 ±9 :E )(4-^)-i w >o. 

Proof. In view of Lemma T2.91 we have w(t,x) > 0, V(i, 2) £ R+ xl. A simple 
calculation shows 

(4-^)- 1 w = -/ e- 2 ^w(t,0c% 



and 



1 r x ^ r°° 

- e~ 2 * / e 2 ^w(t, + je 2 * e~^w(t, 



1 1 f°° 

d x ((4 - d 2 )" 1 ^ = -- e- 2 * J e 2 ^w(t, 0d£ + - e 2 * J e" 2 ^, 
Combining above two identities, we get 

(k 2 ±d x )(4-d 2 x )- 1 w = l (k 2T 2)e- 2 * f e 2 tw(t,0dt 



1 f°° 
+ - (fc 2 ± 2)e 2x / e~ 2 «w(i,£K > 0. 



□ 
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3 Proof of stability 

In this primary section of the paper, we prove the stability theorem (Theorem 
[]} stated in the introduction. Note that the assumptions on the initial profiles 
guarantee the existence of an unique global solution of equation (|1.5[) . The 
stability theorem provides a quantitative estimate of how closely the wave must 
approximate the peakon initially in order to be close enough to some translate of 
the peakon at any later time. That translate must be located at a point where 
the wave is tallest. The proof of Theorem [1] is based on a series of lemmas 
including some in the previous section. 

We take the wave speed c = 1 and the case of general c follows by scaling 
the estimates. 

Note that (p(x) — e~\ x \ € _ff 1 (R) has the peak at x = and 

f 00 2 
E 3 (ip) = / e-^dx = -. (3.1) 



— oo 



Define v u = (4 - d^u = ±e~ 2|2:| * u. Then 

v v {x) = \ f e- 2 ^- xl e~ M dr) = le^ xl -le' 2lxl , Va; E R, (3.2) 
4 Jr 3 6 

and thus 

maxu v = v v (0) = -. (3.3) 

Note ip — d 2 ip = 26. Here, 6 denotes the Dirac distribution. For simplicity, 
we abuse notation by writing integrals instead of the H^ 1 / H 1 duality pairing. 
Hence we have 



E*(<P) = IMIx = / (1 - %)<P (4 - dlY\ dx 

1 



(3.4) 



= 2 / 8(x)(4 - fig)~V(a!) dx = 2v v (0) = ~. 
Jr 3 

Lemma 3.1 For any u S L 2 (R) and £ & R, we have 

E 2 (u) - E2{<p) = \\u-cp(.- 0113c + 4 K(0 " %(0)) > 
w/iere w u = (4 — dl)~ 1 u. 

Proof. This can be done by a simple calculation. To see this, we have 

h - ¥>(• - 0113c = Hx + IM& - 2 /(l - d 2 x Mx - C)(4 - 



Nll + IMIi-4 / <5(x-£)(4-^)-lu0r)cfa: 

Hx + IMIx - 4 ««(0 = IMIx - IMIx + 2\\<p\\x - 
£; 2 (u)- J B 2 (^)+4( % (0)-t;„(0). 
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where use has been made of integration by parts and the fact that E^ip) = 
\\ip\\x = 2u ¥> (0). This completes the proof of the lemma. □ 

In the next two lemmas, we establish two formulas related the critical values 
of v u to the two invariants E%(u) and E%(u). Consider a function ^ u 6 
H s ,s > 3/2 and u > 0. Then < v u = (A-dfj^u G H s + 2 C C 2 . Since 
v u is positive and decays at infinity, it must have n points {£t} i=1 with local 
maximal values and n — 1 points {T]i}™ = i with local minimal values for some 
integer n > 1. We arrange these critical points in their order by 

-oo < £i < ryi < £2 < rj2 < ■■■ < < < & < ^ < ■•■ < »7n-i < < 
Let 

= Mi, 1 < i < n and = rrii, 1 < i < n — 1. (3-5) 

Here, we assume n < +00, that is, there are a finite number of minima and 
maxima of v u . In the case when there are infinitely many maxima and minima, 
the proofs below can be modified simply by changing the finite sums to infinite 
sums. 

Lemma 3.2 Let / u 6 H s 1 s > 3/2 and u > 0. By the above notations, 
define the function g by 

9(X) -\2v u + d 2 x v u + 3d x v u , ZiKxKru, 1 " 1 " n ' (3 ' 6) 

with f]o = — 00 and r] n = -\-oo. Then we have 

P / n n—X \ 

/ g 2 (x)dx = E 2 (u) - 12 [J2 M ? ~ E m * • ( 37 ) 

® \i=l 4=1 / 

Proof. To simplify notations, we use v for v u below. Then u = Av — d 2 v. First, 
we note that 



E 9 (u) 



— / [(1 — <9 2 )u] vdx = / (uv + d x ud x v) dx 
= [ {(Av - d 2 x v) v + [(4 - d 2 x ) d x v] d x v) dx 



Av 2 + 5 (d x v) 2 + (d 2 x 



dx. 



To show p.7p . we evaluate the integral of g 2 on each interval [rji-i, 77J, 1 < i < n. 
We have 

g 2 (x)dx= I (2v + d xx v — 3d x v) dx + (2v + d xx v + 3d x v) dx 
= 1 + 11. 
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To estimate the first term, by integration by parts, we have 

(4v 2 + (d xx v) 2 + 9 (d x v) 2 + Avd xx v - I2vd x v - 6d xx vd x v 

(Av 2 + (d xx v) 2 + 5 (d x v) 2 ) dx-6(v (&) 2 - v (m~if) , 



Vi-l 



where use has been made of the fact that d x v = d x v (rji-i) — 0. Similarly 
II = 



* (Av 2 + (d 2 x v) 2 + 5 (d x vf) dx + 6(v ( m ) 2 - v (&) 2 ) ■ 
g 2 {x)dx = {' (±v 2 + (d 2 x v) 2 + 5 (d x vf^ dx-12v &) 2 +6v {m-if+Gv (^f 



So 

rVi rV- 



'17i-l ''Vi-l 

and 



/ g 2 {x)dx = [ Uv 2 + {d 2 x vf + 5 (d x vf) dx-J2 (l2v (&) 2 - 6v fa-i) 2 - 6v ( m f 



n-1 \ 



= 2fr(«)-12 



m 2 



i=l i=l / 



where use has been made of the fact that v (rjo) = v (rj n ) = and the notations 
in JS3. □ 



Lemma 3.3 With the same assumptions and notations in Lemma \3. 6 A Define 
the function h by 

if \ / -d 2 v u - 6d x v u + 16v u , r)i-i<x<£i, , , . , ^ 

mx) = < i cq i ic ^ 1 < i < n. (3.8) 

wtt 770 = — oo and rj n = +oo. Then we have 

~ / n n— 1 \ 

/ ^(^(ajdi = £7 3 («) - 144 (J2 Mf - . (3.9) 

Jr Vi=i <=i / 

Proof. We still use v for u u . First, note that 



E 3 (u)= / (4v - d xx vY dx = / - (d xx vY + 12 w - 48v 2 d xx v + 64v 3 



dx. 



To show p. 91) . we evaluate the integral of h{x)g 2 (x) on each interval [r7»_i, 
1 < i < n. We have 

h(x)g 2 (x)dx = / (— d xx v — 6d x v + 16v) (2v + d xx v — 3d x v) dx 

i-l Jvi-l 

+ / {—d xx v + 6d x v + 16v) (2v + d xx v + 3d x v) dx 
= 1 + 11. 
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It is found that the first term 



I 1 {- {d xx vf 

Jr/i-i 



+ 12 (d xx v) v + 27d xx v (d x v) - 108vd xx vd x v + 60v d xx v 



54 (d x v) 3 + 216 (d x v) 2 v - 216v 2 d x v + 6Av 3 }dx 



= [ {- (d xx vf + 12 (d xx vf v + 54 (d x v) 3 + 60v 2 d xx v - 54 (O x vf 

Jrii-l 

- I08v 2 d xx v + 6Av 3 }dx - 72 (v (&) 3 - v (^_i) 3 ) 

= - {d xx vf + 12 (d xx v) 2 v - 48v 2 d xx v + 6Av 3 } dx - 72 (v (£) 3 - v 

where use has been made of the following integral identities due to integration 
by parts and d x v (&) = d x v (rii-i) = 0, 

/ d xx v (d x v) 2 dx = ]- [ d x ({d x v f\ dx = 
[ vd xx vd x vdx=f vd x (\- {d x vf\ dx = ~ f (d x vf dx, 
J (d x v) 2 v dx = J d x vd x Q y2 ^ dx = ~\ J v 2 d xx vdx, 
j v 2 d x vdx=[ \d x (v 3 ) dx = i (v(df - v (m-i f)- 

Similarly, 

// = jT' [- (d xx vf + 12 (d xx vf v - 48v 2 d xx v + 64i; 3 ] dx+72 (v ( Vl ) 3 - v (&) 3 ) 



and thus 



f' h(x)g 2 {x)dx= ^ \- (d xx v) 3 + 12 {d xx vf v - A8v 2 d xx v + 6Av 3 

- UAv (&) 3 + 72 (v (t^) 3 + v (^) 3 ) • 
By adding up the above integral from 1 to n, we get 

f h(x)g 2 (x)dx= [ f- (d xx vf + 12 (d xx vf v ~ 48v 2 d xx v + 64v 3 

n n 

- 144$> te) 3 + 72^ (y (v^f + v ( Vl f) 

2 = 1 i=l 

/ n n—1 \ 

= E 3 (u)-lU [Y^Mf-Y,™ 3 . 



dx 



dx 



\i=l 
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□ 

Without changing the integral identities (|3 . T|) and (13. 9|) , we can rearrange 
Mi and rrii in the order: 

Mi > M 2 ■ • • > M„ > 0, mi > • ■ ■ > m„_i > 0. 

Moreover, since each local minimum is less than the neighboring local maximum, 
we have Mi > TOj_i (2 <i < n). The following two elementary inequalities are 
needed in the later proofs. 

Lemma 3.4 For any n>2, assume {Mi}™ =1 and {tth}"~^ are 2n — 1 numbers 
satisfy 

Mi > M 2 ■ ■ ■ > M n > 0, mi > • • • > m„_i > 

and Mi > mj-i (2 < i < n). Then 
(i) 

(Mf - mf_i) < |Mi £ (Mf - mU) ■ (3.10) 

i=2 i=2 

(ii) 

1 1 

(n \ 2 / n \ 3 

M 1 2 +^(Mf-mt 1 )J > U 3 +^(Mf-mti)j . (3.11) 

Proof, (i) For any 2 < i < n, we have 
(Mf - mU) - (Ml - ml,) 

= -- (Mi - rrii-t) (3MiTOi_i + 3MiMj - 2M i m i _i - 2Mf - 2m 2 _ 1 ) < 

since Mi > Mj > mj_i. This implies that desired result in (|3.10|) . 
(ii) Denote 

i i 

(n \ 2 / n \ 3 

Mi 2 + £ ( M ? ~ m li) J and B » = I Mi + ^ (Mf - mf_i) j . 

(3.12) 

We want to show A n > B n (n > 2) by induction. For the case of n = 2, it is 
equivalent to show that 

(Mf + Mf - mf) 3 - (Mf + Mf - to?) 2 > 
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if Mi > Mi > mi > 0. We have 

(Mi 2 + Mf - m 2 ) 3 - (Mf + M| - mlf 

= ZM^Ml - 3Mi 4 m 2 - 2M 3 M 3 + 2M 3 m 3 + 3M 2 M 4 

- 6M 2 M 2 2 m 2 + 3M 2 m 4 - 3M 2 4 m 2 + 2Mf m? + 3M 2 2 m 4 - 2m? 

= (M 2 - mi) (Mi - mi) {3MiM| + 3M?M 2 - 2M x m\ + SMfrrn 

- 2M 2 m\ + 3M 2 3 mi - 2m 4 - 2M? M| + M 2 m 2 + M 2 2 m 2 

- 2MiM 2 m 2 + MiMfmi + M 2 M 2 mi}, 

which is obviously nonnegative by the assumption Mi > M 2 > mi > 0. Assume 
the inequality A n > B n is true for n < k (k > 2). Our goal is to deduce 
Ak+i > Bk+i- Since Ak > Mi > Mk+i > m^, we have 



A 6 - 



fc+1 n 3 



M 2 + ]T(M 2 



(A 2 + M 2 +1 - ml) > (A\ + Ml +1 - ml) (by the n = 2 inequality) 
fc + Ml +1 



> (bI + M^ +1 — ml) (since Ak > -B& by the induction assumption 



= -Bfc+i- 

Thus ^4fc+i > -Bfc+i and A n > B n is true for any n > 2. □ 

The following lemma is crucial in the proof of stability of the peakons. 

Lemma 3.5 Assume uq £ H s ,s > 3/2 and yo > 0. Let Mi = v u (t,^i) = 
max x( zK{v(t, x)} Then for t > 0, 

E 3 (u) - 144B 3 < 18Mi (E 2 (u) - 12A 2 n ) (3.13) 

where u is the global solution of equation il.5\) with initial value uq, v u = (4 — 
<9 2 ) -1 m, and A n and B n are defined in H3.12\) . 

Proof. First, by Lemma 12.81 the global solution u of equation (| 1 . 5[) satisfies 
u(t, x) >0 and y(t, x) — u — d x u > for all (t, x) e R+ X E. We now claim that 
h < 18v for (t, x) G M. + x R. To see this, we rewrite the expression of h as 

hM - J _ + 3dx + 12 ) v - 3d - v + 18v > ni-i <x<&, i < < „ 

«W ~ \ - (a 2 - 35, + 2) v + 3d x v + I8v, ti<x<TH, ~ ~ 

If rji-i < x < £j, 1 < i < n, then > 0. On the other hand, it follows from 
Lemma 12.101 that 

- {d 2 x + '3d x + 2)v = -(2 + 8 X ){A - d 2 x )-\l + d x )u < 0. 

Hence 

- (<9 2 + 3d x + 2) v - 3d x v + 18v < 18v. (3.14) 
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A similar argument also shows that for £j < x < rji, 1 < i < n, d x v < and 

- (dl - 3d x + 2) v + 3d x v + 18« (3.15) 
= -(2 - d x )(4 - 9*) _x (l - d x )u + id x v + 18v < 18v. 

The combination of (|3.14[) and (|3.15|) yields 

h < 18v < 18maxt> = 18 Mi, V (t, x) e K + x R. 

By the notation in (|3.12[) . the integral identities (|3.7p and (|3.9p become 

/ g 2 (x)dx = E 2 (u) - 12A 2 n 

and 

/ h{x)g 2 (x)dx = E 3 (u) - 1MB*. 
Note that when n = 1, A\ = Bi = M\. Relating the above integrals, we get 
B 3 (u) - 1445^ < 18Mi(E 2 (u) - 12A 2 n ). 

□ 

Lemma 3.6 Assume u E H s , s > 3/2 and y > 0. Let 

Mi = v«(i,£i) = max{v(t, x)} 



and 



rn 



i-l 



i=2 



where Mi and wii-i (2 < i < n) are /ocaZ maxima and minima ofv u . If\E2(u) 
E 2 (cp)\ < S and \E 3 (u) - E 3 (tp)\ < S with < 5 < 1, i/ien 

Mi ' 



< VS, 



recalling that v v (0) = h = maxv v , 
(ii) 



< Vs. 



(3.16) 



(3.17) 



and 



(Hi) 



Ws. 



i = 2 



(3.18) 
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Proof. To obtain (i), we first claim that 

1 _ , . . . I 



-E 2 (u)M 1 + — E 3 (u) < 



(3.19) 



In the case when Mi is the only local maximum of v u , we have n — 1, A\ — 
B 1 = Mi and (|3TT5|) follows directly from 1(533} . When n > 2, in view of (|3"TUj) 
and inequality (|3 . 1 0(1 in Lemma [3~4l (i). there appears the relation 

M?-]e 2 (u)M 1 + ±-E 3 (u) 



72 



<lAA[Y J {Mt-mU)-U 1 Y d {Ml 



< 0. 



\i=2 



Define the cubic polynomial P by 



1 



P(y) = y<--E 2 (u)y+-E 3 (u). 



(3.20) 



For the peakon solution, E 2 (ip) = \ and E 3 (tp) = |, the above polynomial 
becomes 



1 



1 



Po(y) = y -7^+w = \y-a) \y+o 



12" 108 



1 



(i 



1 



(3.21) 



Since 



Po(M x ) = P(Mi) + i (Ba(u) - £ 2 (^)) Mt - ± (£ 3 (u) - E 3 (<p)) 
and P(Mi) < by lETTSj) , it follows that 



Mi 



G 



< I { e 2 (u) _ E 2 {ip)) Mi - 1 (£ 3 (u) - £ 3 M) 



On the other hand, observing E 2 {u) — 12A^ > 0, we have 

< Mi <A n < v/£ 2 (u)/12 < ^(1/3 + (5) /12 < i 
when <5 < 1. It is then inferred from (|3.22[) that 
1 



(3.22) 



(3.23) 



Mi - 



< yJl\E 2 (u)-E 2 (i P )\ + ^\E 3 (u)-E 3 ( ( p)\<V5. 



We now prove claim (ii). When n = 1, Ai = Mi and it is reduced to (i). 
When n > 2, it is thereby inferred from (|3.13p that 



_ -^(u)^ + — £ 3 (n)<0, 



(3.24) 



due to < Mi < A n and < B n < A n by Lemma [3.41 (ii). In consequence, 
(|3.17p follows from the same argument as in part (i) . 
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(iii) can be obtained from (i) and (ii). In fact, combining (i) and (ii), we 
have 

~ A n + Ah 

which implies (|3.18p by using (|3.23p . □ 

Proof of Theorem ED Let u e C([0, oo), H s ), s > 3/2 be the solution of (|L~5|) 
with initial data u(0) = uq. Since E 2 and E3 are both conserved by the evolution 
equation (|1.5j) . we have 



E 2 (u(t, ■)) = E 2 (u ) and E 3 (u(t, •) = E 3 (u ), Vt > 0. (3.25) 
Since ||«q — <p\\x < e, we obtain 

\E 2 (u ) - E 2 {y)\ = \{\\u \\ x - |M|x)(||u ||x + |M|x)| 
<£(2\\<p\\ x +e) = e(^=+e^ <2e, 

under the assumption e < \ . In view of (]3 . 25[) , the assumptions of Lemma 13.61 
are satisfied for u(t, •) and 8 — 2e. It is then inferred that 



< V2e, Vt > 0. (3.26) 



By (|3 . 25[) and Lemma f3.ll we have 
\\u(t, •) - <p{- - Zi{t))\\ x = E 2 (u ) - E 2 (<p) + 4(^(0) - ««,(*, (*))), Vt > 0. 
Combining the above estimates yields 

|K*, •) - (p(- - < V2£^4V27 < 3e 1/4 . Vt > 0. 

Estimates (|1.9p and (jl.lOp then follow directly from Lemma [3.61 (ii) and (iii). 
This completes the proof of Theorem [TJ □ 

Remark 2 We make several comments. 

(1) By US. 23\) . Mi = maxn„ < \J E 2 (u) j\2. For peakons ctp, we have 
max.v clp = \J E 2 {op) j\2 = ^c. So among all waves of a fixed energy E 2 , the 
peakon is tallest in terms ofv u . 

(2) In our proof, we use inequality H3.22\ ) to get estimates US. 16}) and \3.18} ) 
more directly, compared with the argument in [12] by analyzing the root struc- 
ture of the polynomial P{y). Moreover, it implies that the peakons are energy 
minimizers with a fixed invariant E3, which explains their stability. Indeed, if 
^3(11) = Es^ip), it follows from i3.22\) that E 2 (u) > E 2 ((p). The same re- 
mark also applies to the CH equation and shows that the CH-peakons are energy 
minima with fixed F3 . 

(3) Compared with \12]j , our construction of the integral relations {3. 7f ) and 
is more delicate. It not only is required in our current case, but also 
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provides us more information about stability via M.10\) . For the CH equation, 
even if the orbital stability is proved by a simpler construction \12f . our approach 
can also give the additional stability information. More specifically, for the CH 
equation il.3\) with yo > 0, by refining the integrals of '12, Lemma 2] to each 
monotonic interval of u, one can obtain 

F 3 (u) - ~Bl < M 1 (F 2 («) - 2A 2 n ) , 

where F% and F 3 are defined in \1. 7| ), and A n and B n in h3.12}) with Mi and rrii 
being the maxima and minima of u, respectively. Hence, estimate \1.10\) may 
be obtained by following the proof of Lemma \3.6\ 
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